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Abstract 

In one of the very few exact quantum mechanical calculations of fugacity coefficients, Dodd and 
Gibbs (J. Math.Phys.,15, 41 (1974)) obtained 62 and 63 for a one dimensional Bose gas, subject to 
repulsive delta-function interactions, by direct integration of the wave functions. For 62, we have 
shown (Mol. Phys., 103, 1301 (2005)) that Dodd and Gibbs' result can be obtained from a phase 
shift formalism, if one also includes the contribution of oscillating terms, usually contributing only 
in 1 dimension. Now, we develop an exact expression for 63 — 63 (where 63 is the free particle 
fugacity coefficient) in terms of sums and differences of 3-body eigenphase shifts. Further, we 
show that if we obtain these eigenphase shifts in a distorted-Born approximation, then, to first 
order, we reproduce the leading low temperature behaviour, obtained from an expansion of the 
two-fold integral of Dodd and Gibbs. The contributions of the oscillating terms cancel. 

The formalism that we propose is not limited to one dimension, but seeks to provide a general 
method to obtain virial coefficients, fugacity coefficients, in terms of asymptotic quantities. The 
exact one dimensional results allow us to confirm the validity of our approach in this domain. 

Keywords Virial Coefficients • Hyperspherical basis ■ Boson systems 
PACS 05.30.Jp 05.70.ce 34.50.S 
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Introduction 



As demonstrated in NaturepQ the behaviour of quantum models, even in one dimension, and of 
such gases, remains of current interest! In that paper, the paramount interest is the experimental 
investigation of a system that is not ergotic and does not equilibrate after thousands of collisions. 
Still, the description of a quantum mechanical gas, in equilibrium, needs adequate description. So, 
for example, we note the fairly recent article by Hussien and Yahya [2], on the quantum corrections 
to the classical fourth virial coefficient for square well potentials. 

Our efforts are, instead, directed to the continuing development of a fully quantum mechanical 
formalism, applicable to the higher virials, and especially useful at low temperatures. 

Let us be clear, there is an enormous difficulty in evaluating the required traces in Statistical 
Mechanics, which involve 3 or more particles. Our adiabatic hyperspherical formalism, uniquely, 
can 'deal' with 3 or more quantum mechanical particles, whether for Boltzmann or Bose or Fermi 
statistics, in arbitrary dimensions. This is very important: it can be extended to more particles. The 
hyperspherical formalism can, in principle, be used to represent an arbitrary number of particles. 

Another 'key' element is that we express our results in terms of 'asymptotic' quantities, in our 
particular case here: 'eigenphase shifts'. This can be done for the 2-body problem, and the essential 
'jump' is from two to three and more particles! 

We are continuing to extend our formalism, and we use 1-dim to validate our results and to verify 
that we can carry out our fugacity program. In this paper, we show how we 'deal' with Bose statistics 
and repulsive forces. Previously, we have limited ourselves to Boltzmann statistics, but developed the 
present approach which relies on the use of a complete hyperspherical adiabatic basis (3J. 
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We recall that, many years ago, Dodd and Gibbs [I] were able to evaluate the second and third 
virials of identical Bose particles in one dimension, subject to repulsive delta function interactions. 
To do this, they introduced a complete set of energy eigenfunctions in the traces, and then integrated 
to obtain answers in terms of single and two-dimensional integrals. 

Their original aim, as is clear from their paper and from Gibbs' thesis, was to obtain a formalism 
in terms of the S-matrix, following the formalism of Dashen, Ma and Bernstein [5 , in an effort to 
obtain the virials in terms of scattering quantities. This was not possible, as the authors determined 
that the "formal limiting processes are not valid for the singular amplitudes of this system". 

An alternative formalism, proposed by Larsen and Mascheroni [5J, was a step aimed at ultimately 
expressing the higher virials in terms of eigenphase shifts and bound state energies, and, if necessary, 
other on-shell properties, in a hyperspherical harmonic formalism. Under some constraints, one was 
then able to recover the classical results for the 3rd virial, through a semi-classical approximation 
[TJ. Low temperature calculations have been carried out for repulsive potentials in 2 dimensions [8 , 
and a more sophisticated formalism, potentially able to 'handle' correctly the bound states, has been 
presented for Boltzmann statistics [3 . It involves the use, from the beginning, of a hyperspherical 
adiabatic basis. Here, we extend that formalism to apply to Bose gases. 

To summarize, our aim, in this paper, is to use the excellent results of Dodd and Gibbs to confirm 
our virial formalism for the third fugacity coefficient, for this particular dimension and interaction, 
for which we have analytical and numerical results [SJ [TUl [TT] . More specifically, we will show that the 
relevant 63 — 63 can be rewritten as an integral, over the energy, of a Boltzmann factor together with 
the sums and differences of eigenphase shifts, for 3-body wave functions in one dimension. Using, then, 
a distorted-Born approximation for the eigenphase shifts, we reproduce the leading low temperature 
behaviour, obtained from an expansion of the two-fold integral of Dodd and Gibbs. 
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Statistical Mechanics 

In this section we start with the Grand Partition function in 3 dimensions: 



Q = J2 z " Tr (e"^") = 1 + zTr (e^ 1 ) + z 2 Tr (e^ 2 ) + z 3 Tr (e^" 3 ) + 



(1) 



n=0 



The fugacity z equals exp (/i/reT); j3 = 1/kT, where fi, k and T are the Gibbs' function per particle, 
Boltzmann's constant and the temperature, respectively. H n and T n are the n-particle Hamiltonian 
and kinetic energy operators. 

We note that no factorials (1 /factorials) appear in this development. This is correct for Bose or 
Fermi statistics. Also important, is the fact that 

Tr (e _/3H ") — > V n as the volume V — > large, 

leading to the divergence of the individual traces in the thermodynamic limit. If we take, however, 
the logarithm of the Grand Partition function Q, we obtain 



InQ 



zTr (e^ Tl ) 



+z' 



Tr (e-^)-i[Tn 



Tr (e-l 3 ^) ~ Tr (e^) Tr (e"^) + 1 [ Tr ( e "^)] : 



(2) 



which, when divided by V, gives coefficients of the powers of z, which are independent of the volume, 
when the latter becomes large. They are the fugacity coefficients bi . We can then write for the pressure 

and the density 



p/ K T=(l/V) \nQ = J2 b 



1 

N/V = p = J2 lb i z '> 
1 

and the fugacity can be eliminated to yield the pressure in terms of the density, in the virial expansion: 

00 



One dimension, three Bose particles 



We follow the procedures outlined in the 'Bogoliubov' lecture[3], developed for Boltzmann particles, 
and adapt the formalism to Bose Statistics. 

From Eq.([2]), and changing from the volume V to a length L, an expression for the 3rd fugacity 
coefficient can be written as 



Tr (e 



Tr (e- m )Tr(e-^) +-Trl 

3 



-/3Ti 



(3) 



Note that each of the 3 terms in the bracket, above, grows as L 3 asymptotically, for large L. Therefore, 
subtractions must take place, so that 63 might become L-independent, as is required. If we factor out 
the contribution of the CM. (proportional to L), each term will have a dependence proportional to 
L 2 . 

To help us decrease the L dependence, let us subtract, from each of terms above, the equivalent 
term without interaction. 

Subtracting, therefore, 63 from the 63, we obtain: 

63 - b° 3 = j [Tr (e-? 11 * - e-^) - Tr (e~^) Tr {e~? H > - e~^)} (4) 
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Adiabatic Preliminaries 



For the 3 particles of equal masses, in one dimension, we first introduce center of mass and Jacobi 
coordinates. We define 



V 



1/2 



(xi - x 2 ) 



2\ 1/<2 / xi + x 2 



X3 



R = g (Xl + X 2 + X 3 ) 



where, of course, the Xi give us the locations of the 3 particles. This is a canonical transformation 
and insures that in the kinetic energy there are no cross terms. 

The variables £ and r\ are involved separately in the Laplacians and we may consider them as 



acting in different spaces. We introduce a higher dimensional vector p 



and express it, here 



in two dimensions, in a simple polar coordinate system (using the radius p and the angle We write 

f] = pcosi?; £ = psiwd, —it < i9 < 7r; < p < oo. (5) 

If we factor a term of p 1 ! 2 from the solution of the relative Schrodinger equation, i.e. we let 
tp — 4>/ p 1 / 2 , we are led to: 

d 2 



where 



Hp — J 

P 2 



dd 2 



2m 



m is the mass of each particle, k = y/2mE /K 2 and E is the relative energy in the center of mass. 
We now introduce the adiabatic basis, which consists of the eigenfunctions of part of the Hamiltonian: 
the angular part of the kinetic energy and the potential. 



H p B t {p,0) = At(p)B e (p,0), 



(6) 



where t enumerates the solutions. 

Using this adiabatic basis, we can now rewrite the Schrodinger equation as a system of coupled 
ordinary differential equations. For a typical solution, we write 



(7) 



and obtain the set of coupled equations 
d 2 



_ -A e (p) + k 2 ) MKp) + 2^Cf(p)—Mk,p) 



(8) 



where we defined: 



Cf(p) 

dUp) 



I d&BW,p)^-B t ,W,p) 

/rj2 
d0BW,p)-^B t ,(0,p). 



(9) 
(10) 
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Three binary interactions: case (123) 



In our model calculation, the interaction between our particles will be a repulsive delta-function, with 
'strength' g > 0. For 3 particles, we will thus have an interaction potential V{p, §): 

V(p, tf) = g[S(V2p\ costfQ + 6(V2p\ cos(i? - 2tt/3)|) + S(V2p\ cos(tf + 2tt/3)|)]. (11) 

In a set of previous papers [9-15] we, and colleagues, have studied these systems (with g positive 
and negative), calculated eigenpotentials, adiabatic bases, phase shifts, eigenfunctions and binding 
energies - where appropriate. We refer to the Appendices for a selection of results, that we need in 
our present calculations. 

The adiabatic eigenfunctions Bk{p,$) can be characterized by their symmetry, together with the 
index K, as they reduce to simple normalized surface harmonics when the potential V(p, 1?) — >■ and 
their associated eigenvalues Ak — > (K 2 — j)/p 2 - 

For these surface harmonics, in the completely symmetric 1-dimensional representation Ti, we 
need [12] to set K = mod 3 . The harmonics will then be: 

for K — even, -4= cos(K-d) with —}= for K = 

for K = odd, 4= sm(K-d) y ' 

[The completely antisymmetric 1-dimensional representation T\ requires, in the above, to inter- 
change the role of the sines and cosines (clearly without the constant term). Finally, in cases where 
K 7^ OmodS, the sines and cosines form pairs belonging to a 2-dimensional representation T2. ] 

The Bk (p, can then be constructed from the harmonics or obtained as solutions of a differential 
equation, Eq.©. The resulting forms are then|14|: 



for K = even, A K (p) cos(q K (tf - *f )) where q K satisfies q K tan(^) - ^ 



for K = odd, A K (p) sm(q K (i9 - ^)) where q K satisfies q K cot(^) ^ -E££ 



(13) 



7T 

_ 6 

and where m is an integer such that \d — mn/3\ < tt/6 and c = (2m/H 2 )(3g /tt^/2). 

For each value of p, we need to solve the transcendental equation for qx- We can do so numerically, 

or alternatively seek expansions in p or 1/p, for small or large values of p, respectively. 

Ak{p) is simply the normalizing factor, such that the Bx{p, $) are orthonormal when integrating over 

the angle, with a choice of integrating over a 'sector', say — 7r/6 to 7r/6, or over a more extended range 

such as from to 2it. 



Two interacting particles + a spectator: case (12,3) 

In this case, which we also require, we use our 3-body formalism to treat three particles, only two of 
which are interacting, the other one playing the role of a spectator. It is very important that this be 
done in a way closely paralleling the treatment for the three interacting particles, because we will need 
to subtract traces, associated with lengths (volumes in 3-dimensions) and they must be evaluated in 
a similar way in all of the terms of our expressions. 

Since our need is the evaluation of traces, we can choose particles 1 and 2 as interacting, particle 3 
as the spectator. Letting -d — > — -i9 + 7r, and therefore letting 77 go to —77 and £ go to £, interchanges the 
role of particles 1 and 2. The symmetrical combinations for the harmonics (i.e. which treat particles 
1 and 2 in a symmetrical fashion) will then be: 

for K = even, A= cos(Kd) with— h= for K = 

for K = odd, -±= sm(K<&) { ' 

where, of course, there is no longer a restriction of mod 3. K now runs over 0, 2, 4, 6... or 1, 3, 5.... 
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The B^ 3 (p,d) will read: 



for A = even, A 1 ^ 3 (p) cos(g r if(i? — mir)) where qx satisfies qx ta.n(^^) = — 

(15) 

for A = odd, A 1 ^ ,3 {p) sm(qK (t? — m7r)) where g^- satisfies <7_r- cot(^p) = — 
and where m is an integer such that |i? — to7t| < 7r/2 and, again, c = (2m/H 2 )(3g/n\/2). 

The A]^' 3 (p) are the normalizing factors, such that the B K 2 ' 3 (p,'d) are orthonormal, when inte- 
grated over the angle - with a choice of integrating over a 'sector' (here the right-hand or left-hand 
part of the range) or over the whole range of to 2tt. 

The qx, here, are really q^' 3 and, of course, as in the fully interacting case, are also functions of p. 

A's, C's, D's 

Given the expressions for the 2?'s ... we can then obtain expansions for the behaviour of the qx as 
functions of p, and therefore for Kk and for the C's and the Z?'s .... For us, the most important 
behaviour is that as p becomes large .... i.e. expansions in l/p. 

We then find, [see Appendix B], that for p — > oo the effective centrifugal term (i.e. the l/p 2 term) 
reads : 

Ak = (J ^^- + ■■■ for the case (123) 

(16) 

A 12,3 = (K+lf + for ^ ^ ^ ^ 

and in both cases, the C and D terms behave as : 

f + A^A' 
^ a \ ^ + ... A = A' 

(17) 



d 

U K do 



+ A^A' 
A = A' 



In both cases, to obtain a non-zero result, the indices of the A and A' matrix elements must both 
belong either to an order which is even or to an order which is odd. There are no cross terms. 

Coupled equations and asymptotic behaviour 

Given N coupled equations, we will have N independent solutions (well behaved at the origin) of the 
coupled differential equations. These can be expressed as: 

4> K {k,-fi) = Y i BK>{p,#)<t%.{k,p) (18) 

K' 

where the <p^, - the amplitudes - can be chosen to have as asymptotic behaviour, either: 

4>K>(k,p) = 6%, y/k~p JK>+ 3 (kp) + Wg,(k) ^N K , +3 (kp) (19) 

or _ 

4$,(k,p) = 5%, Jkp J K >+i{kp) + W K ,(k) yfkp N K , +1 {kp), (20) 

as appropriate for the desired case. 

[I.e. starting from amplitudes with A^, and B^, as coefficients of the Bessel and Neumann functions 
(multiplied by \fkp), we take linear combinations of the solutions - multiplying by the matrix A^ 1 - 
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to obtain equations (fT!?]) and (f2"U)) .] The upper index K, chosen to enumerate the solutions (|18|) . is 
associated with the Jx{kp) which dominates the solution for small values of p, and takes on the same 
set of values as the lower index K' . 

The wave functions and the amplitudes can be chosen to be real, and the and the W K will be 

real and symmetrical. They can be diagonalized to yield tangents of eigenphase shifts — tanJ^-(fc) 
and — tan ((5w) K )(k), respectively. 

Further, we can then take linear combinations of these solutions, denoted as <f) v K , , such that for the 
fully interacting system : 

r K >(k,p) -> (M V2 C£'(fc) [cos^(fc) JK> +3 (k P ) - smS^(k) N K , +3 (kp)} (21) 

where K' denotes the K'th amplitude of the solution v and the eigenphase shifts are associated with 
the matrix W. In the fully interacting case, these phase shifts will only differ by — 3tt/2 from those of 
the genuine iS-matrix for this problem. (I.e. — 3ir/2 + Sw = Ss) 
In the (12,3) case, we obtain similar expressions. 

We note that since the L K' = even' amplitudes couple only with the L K' = even' amplitudes 
and, analogously, the l K' = odd' amplitudes couple only with the l K' = odd' amplitudes, we can 
therefore separate our set of solutions 4> V K , into 2 sets, one with amplitudes involving only even values 
of K' and the other with amplitudes involving only odd values of K' . 
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The phase shift expressions for the Traces 



1st term - with the totally symmetric adiabatic (123) basis 
Consider the first term of the Eq. (|4]) : 



31/2 
At 



Tr 1 



(e-^ ei ) -Tr(e-^ e! )l , (22) 



seen here, after we have separated the Hamiltonian into two parts, one corresponding to the center 
of mass coordinate and the other to the relative motion, and taken the trace of the center of mass 
contribution. The latter contributes a factor that depends on the temperature At = (/i 2 /27ttokT) 1 / 2 
and a length that cancels the factor 1/L shown in Eq.(U). The first term of the above can now be 
written as 



31/2 
At 



/poo 
d-f J dk (k,^)^ {k,jf)e- f! ( h2/2m ) k2 (23) 

where we introduced a complete set of energy eigenfunctions, defined over the hyperspherical space '. 



As we have indicated before, we let ip u — (jf / p 1 / 2 . and then expand the eigenfunction <j> v in terms 
of an adiabatic basis B^? 3 , which is completely symmetric, complete and normalized over each of the 
intervals (m7r/3 — 7r/6, rmr/3 + 7r/6), where m = 0, 1, • • • 5. Each of these intervals is associated with 
one of the 6 possible permutations of the coordinates of the system. Integrating over just one of the 
sectors is to integrate over all of the possible configurations of the physical particles in a Bose system. 

Alternatively, if we wish to integrate over the interval — n to 7r we can divide the angular functions 
by VS\, taking into account their six- fold symmetry, and obtain a symmetric adiabatic base Bk (p, 0), 
normalized over this extended angular interval. Expanding the functions in terms of the symmetric 
adiabatic base, expression ([23)) now reads: 

^j™dkj °°dp £ <p%(k,p)r K ,(k,p) 

v,K,K> (24) 

j% d$B* K (p, 0) B K , (p, 0) e -^ h2 l 2m > 2 

Here v stands for the i/th solution and K, which takes on only even or only odd values for a given v, 
labels the components. Integration over then gives 



gl/2 i-oo 

At Jo 



p OO p oo 

/ dfc]T / dp^ K *(k,p)^ K (k,p)e- f3 ( h2 / 2m ) k2 (25) 
Jo .. ^ Jo 



We note that had we expanded in terms of B^? 3 (p, 0), and integrated over one of the sectors, we 
would have obtained the same expression. 

The radial functions satisfy the coupled radial equations: 

-^-A K (p) + k 2 ^J 4? K {k,p) + 2^C£'(p)A^,( M 



££>f'(p)«M=0 (26) 



A"' 



and we can choose our solutions to be real. 

To evaluate the integral of the square of the wave function over p we now use a procedure similar 
to that used in our previous papers[l6]. We integrate the square of the wave functions within a sphere 
of radius p ma x, and then let p max go to infinity. We use the Green's function 'trick' [3], and obtain 
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fc 2 - (fc') 



K 

*&(M f-/K{k',p)~<f> V K{k',p)^ K {k,p)] 



K 



dp 



evaluated at p = p.„ 



I.e. our identity is: 



K 

+ {P-(kr)^K(k,p)r K (k',p) 

K 

+ 2 E i p) °k (p) «v(M 



(27) 



(28) 



and we integrate with respect to p. Using then the fact that <jf K goes to zero, as p itself goes to zero, 
and that Cj£ decreases fast enough for p large, we are left with the expression displayed earlier (that 
of our 'trick'). 



We now put in the real version of the asymptotic form of our solutions, Eq. ([2"Tj) . oscillatory solutions 
valid for p max large, and use l'Hospital's rule to take the limit as fc' — > k. Thus for the interacting 
system 

<j> v K (k,p)^(kp) 1 / 2 C» K (k)[cos6^(k) J K+3 (kp)-Bm6k(k)N K+3 (kp)]. (29) 

Here, we note again that we come to a major difference with what happens in other dimensions. We 
have seen, Eqs. (|16|) . that for large values of p, the effective asymptotic 'centrifugal' term involves 
K + 3, instead of the K that appears for small values of p. The corresponding 'free' particle solution 
can be chosen to be: 

{4Py K {k,p) = {kp) 1 /*c» K {k)j K {kp) 

We note that 8 V is associated with the solution v, and appears in each of the amplitudes of this 
particular solution. The coefficient C V K , the mixture coefficient, tells us how much of each normal- 
ized amplitude appears in the solution. These two quantities therefore completely characterize the 
asymptotic behaviour of the j/th wave function. We expect that 

E(Q)*Q = i. 

K 

(Or the real equivalent.) 

Inserting this into our integrals, and taking the limit kl — > fc, we find that 

l^xik, P)\ 2 dp -> --jT^ik) + -Pmax + osc. terms 
and, thus, that 

E l PmaX (\<t>K(k, P)\ 2 - \{<f)K{k, P )\ 2 ) dp -> + osc. terms (30) 

We let pmax go to infinity, and usually the oscillating terms do not contribute to the integral over the 
energy. In 1 dimension, however, we have to be more careful. 
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Using the asymptotic form of real amplitudes 



-C£(jfe)sin 



hp- ~(K + 3)?r + -TT + (A) 



(31) 



and l'Hospital's rule, when taking the limit as k 1 
can be written as 



k in Eq. ([27p . we find that the integral over p , 



S [P™* + (*)] + 5fE ^'k( k f I sin [ 2 ( k Pma X ~ \{K + 3)tt + ±tt + S" (k))] (32) 



When we subtract the term without interaction from the above expression, the integral over k takes 
the form, involving the eigenphases: 



and the contribution from the oscillating terms: 



(33) 



C Jo°° dk ^J2(°K( k )) 2 {* («n [2 {k Pmax - \{K + 3)tt +|7r + ^ (fc))] 

(34) 

3in [2 (kp max - ii^Tr + ±7r)])} e-^VH* 3 



— sm 

which can be rewritten as 

3 1 / 2 poo 



C /o°° dfc E W" { SE ( cos t 2 t cos (2^ (*0) + 1] 

(35) 

- sin [2 (kp max )} sin (2^ (*)))} e -/^ 2 /2™)* 2 



where we have used the fact that the K are either all even or all odd for a given eigenfunction v, 
giving a value of —1 to =p, when the K are even and +1 when they are odd. For a given v 1 we can 
therefore sum over K. In fact, since we have isolated the K dependence of the trigonometric terms, 
we can sum the (C^-) 2 , over K . 

We should also mention that when we perform the subtraction of the 'free' wave function from the 
interacting one ... we must choose precisely the same C's for the free solution, as the ones that are 
mandated when we have an interaction! We are however free to do so! 

The phase shifts, that appear in the expressions above, are the VF-phase shifts ... that go to 
zero as k tends to zero. In our last equation, therefore, the leading term of cos (25) + 1 will give +2 
... and the individual integrands will 'blow' at the origin. We have a major problem! Unless our 
reasoning is incorrect, our only hope is that alternating signs and other similar oscillatory terms - 
arising from the (12,3) contributions - will enable us to precisely cancel what we have here. The term 
sin [2 (kp max )] sin (28 v (k)) jk will not contribute to the /c-integral, when we take the limit p max — > oo 
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2nd term - with the symmetric adiabatic (12,3) basis 

We take the second set of traces in Eq. ^ and rewrite it as an expression involving a trace of a 3-body 
operator: 

Tr (e-^ Tl ) Tr (e^ 2 ) = Tr (V^ 2+Tl >) (36) 

This is an essential trick. To obtain the L-independence of 63, we need to subtract related quantities 
from all of the traces, i.e. evaluate them with same or similar harmonics or adiabatic functions, have 
sums which, though they might individually diverge if summed to infinity, can be subtracted to yield 
finite results. They can each be evaluated for a finite number of terms, subtracted by similar sums, 
and then permitted, for a fixed energy, to converge [This has worked for a Boltzmann calculation^ ] 
or - in our present case - subjected to an accelerated convergence procedure, subtracted and found to 
yield a finite result. 

We now expand the wave functions, eigenstates of a 3-body Hamiltonian of 3 particles, for which, 
say, particles 1 and 2 are interacting, with particle 3 acting as a spectator, in terms of adiabatic 
angular functions and p-dependent amplitudes. 



J2 b k' 3 mTk (m 

K 



(37) 



[We note that v and K are not the same indices as appear in the previous section, and do not assume 
the same sets of values.] 

The angular functions, B 12,3 , will be symmetric upon interchange of particles 1 and 2, but will 
have no symmetry with respect to interchanges between particles 1 and 3 or 2 and 3. They will be 
eigenfunctions of Eq.©, with the potential V equal to: 



V = gS (\xi -x 2 \)= gS (V2p\ cosi?|) 



and normalized over 1/2 of the $ space, for example from — ir/2 to +7r/2, since this corresponds to a 
complete sampling of the possible physical configurations. Alternatively, we can define a Bi equal to 
B 12 ' 3 divided by \/2, which can be integrated from — tt to +tt and still be normalized. The equivalent 
to Eq.([24]) is then 



31/2 
At 



dke-^ 2 ' 2 ™)* 2 



dp 



E 

>,K,K' 



4>k (k, p) <\> K > (k, p) / d-d B K (p, 0) B KI (p, 0) 



(38) 



and the equivalent of Eq. (|25|) is then 

31/2 



ol/z roo poo 

i- / dke^( h2 / 2m ) k2 y dpj* 



(k,p) <j) K (k,p) 



(39) 



We can now subtract a trace with 'free' particles, with wave functions that are only symmetric 
under an interchange of particles 1 and 2, integrate over p and define eigenphase shifts 6 (fc) for this 
3-body system, following steps similar to those of the previous section. There are differences, though. 
Instead of the K + 3, which appeared in the asymptotic behaviour of Kk, we find that the effective 
asymptotic 'centrifugal' term involves K + 1. The phase shifts must therefore be defined in terms of 
Bessel and Neumann functions of order K + 1. They are, again, VF-type phase shifts ... which go to 
zero as k goes to zero. 

For the contribution of the 2nd term, we then arrive at: 



-P{ji 2 /2m)k 2 J. 



(40) 
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where the sum over v must involve the even and odd values of K, but a particular ip will only involve 
the even or the odd values of K. The oscillating terms yield: 



Jo°° dk ^J2 [ C K(k)\ 2 {i (sin [2 (fc/w - + 1 ) 7r + ^ + <T (fc)) 
- sin [2 (fc/w - IJTtt + Itt)])} e-^C 2 /^)^ 
which can be rewritten as 

C /o°° dfc E (T), (cos [2 (fc,w)] [cos ( 2 r (A;)) + 1 
- sin [2 (k Pmax )] sin (2*" (fc))) } e~^' '' /2m > 2 



(41) 



(42) 



where, as before, =p takes the value of —1 when the K are even and +1 when they are odd. 

We see that, just as in the previous oscillating terms for the 1st term of 63, again, the cosine term 

poses problems and the sine term does not contribute to the fugacity coefficient. 
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Distorted Born 



From the foregoing, we see that to obtain analytical results for 63 — 63, we shall need an expansion of 
the eigenphase shifts, for small values of k, our customary wave number. 

The 'conventional' result for a 2-body potential, as stated by Mott and Massey [17] . is that given 
an orbital angular momentum £H, and a radial potential that falls off as r~ s for large r, then 

lim k 2l+1 cot(^) = constant 

if £ < (s — 3)/2 . This is the usual result, in three dimensions, for short range potentials. 
If, however, I > (s — 3)/2, Mott and Masey then obtain 

lim k s ~ 2 cot(np) — constant 

a very different behaviour, that they show can be obtained by a Born approximation. The point is 
that, here, it is the long-tail behaviour of the potential that drives the low-energy behaviour of the 
phase shifts. 

In our hyperspherical formalism we see, for repulsive delta-function interactions, from our equations 
(fTO)l . that when p is large, the role of the I + 1/2 is played by K + 3 or K + 1, respectively and, 
therefore, the criterion for the dominance of the 'tail' becomes (K + 3) > (s — 2) for the (123) case 
and (K + 1) > (s — 2) when discussing the (12, 3) physical problem. We note that this is satisfied for 
our adiabatic potentials. 

Still, the recommended use of a Born approximation is not so straightforward. This is due to a 
peculiarity arising from working in one dimension. In a manner not arising in other dimensions, for 
any of our equations, the effective centrifugal term changes as a function of p, due to the interaction. 
Thus, for the (123) case and a given K, for 'low' values of p it assumes the centrifugal value of 
(K 2 — l/4)/p 2 , while for 'large' values of p, thanks to the interactions, it takes on the asymptotic 
value of (K + 3) 2 /p 2 . 

Thus a 'conventional' Born iteration, with the 'low-p' expression for the centrifugal term, and using 
Jtcikp) as a zeroth order approximation for the solution, would yield iterates that would not form 
an expansion in powers of k, but would in fact each contribute to a constant term in a power series 
expansion of an R- matrix. [Here, to make the point, we just discuss one equation.] 

We were then led to the development of the 'Distorted Born' |T5l [T31 HP] . 

For ample discussion, and details, see J.J. Popiel's thesisQS). Below, we follow his general outline, 
for the case (123), which we have extended to include Ks that are odd, and also to the physics that 
we denote by (12,3). We are happy to have recovered and extended his results; we simply note that 
we differ by a few signs (which, in some crucial equations, compensate!). 

For low energies, we write: 



d 2 (K + 3) 2 



4 



dp 2 p 2 



#(M = £( )°° (43) 

K" 



where 



( )k " (p) = I A* (p) - v " y ) <5£ - Vk (P) 2C£ (P) f ■ (44) 

The superscript 00 on ( )°° is meant to imply that we seek to use the large p expansion of ( )j£ . 
Asymptotically, it will behave as 1/p 3 . The matrix C, which only has off-diagonal elements, is propor- 
tional to 1/p 2 , but only appears in the coupled equations as C j- and hence for large p is effectively 
proportional to k/p 2 . The K' denotes the K'th solution. 

An iteration of our coupled equations will yield a power series in k with each new term in the 
successive iteration beginning with a higher power of k. 

We now write an integral equation, valid for small values of k: 



(k,p) = (<P%')°(k,p) + I G K+3 (p,p';k) 53 ( )°° <f>%'„(k,p')dp' (45) 

K" 
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where [<f>^ )° is the leading term in an expansion of for low energy. Our choice is 

(0f )° (k, P ) = 6%' Jkp J K+ 3(k P ) (46) 

and3 

{ +Jk V^P J K +3(kp)VW N K + 3 {kp') p<p' 

G K+3 (p,p';k)= I _ _ . (47) 

I + 5£ Vkp N K+3 (kp)^ J K+3 (kp') p> P ' 

With this, the large p behaviour of 0^ is 

<t>K (M = ^" Vfcp -/*+ 3 (M + Wg'(k) Vkp~ N K+3 (kp) (48) 

where 1 

W f'W = +2l / VkpJ K+3 (kp') r <t>K»{k,pi) dpi . (49) 

We note that the matrix W now contains all of the scattering information beyond the (k — > 0) result, 
and completely characterizes the asymptotic form of the wave- function.. 
Iterating once, gives us: 

tK' 



Wg (k) = 1st Distorted Bom + ... (50) 



where 1 

1st Distorted Born 



( h 2 k\ ( 24y/2 (K+3)(K'+3) \ 

V mg ) \ tt [(K+K'+6)^-l] [(K-K')^-l] / 



-<5£ x ^ 



. 2 



(¥ + T2) (K+3) 9 1 (51) 

(K+2)(K+4) Att{K+2)(K+3)(K+A) f 



0(k 3 ) 



The 2nd Distorted Born yields a k 2 contribution which exactly cancels that obtained from the 1st 
Born! So, the total contribution of order k 2 is zero. 

For the case (12,3), and W K , the corresponding formula is: 



1st Distorted Born 



h 2 k \ f jvg (K+1)(K' + 1) \ 

mg J \ tt [(K+K'+2)*-l] [(K-K') 2 -l] ] 

(52) 



0(k 2 



The terms of 0(k 2 ) are very much more difficult to evaluate, as the individual diagonal terms for 
K = diverge! As we can see from Mott and Massey's discussion we are, for K = (and s > 3, 
corresponding to higher inverse powers appearing also in the expansion of the interactions at large 
distances) , at the limit of the domain of validity of this method of calculation. 

In Appendix D we show for a simplified version of a model mentioned in K. Chadan and P.C. 
Sabatier's book |19j . for which the exact solution is known, that our Distorted Born approach gives 
exactly the correct coefficients, where the Mott and Massey criterion states that it should (and which 
is also precisely where the integrals of the expansion converge). 



1 Change of sign from Popiel's result 
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The Phase Shift Sums 



from Wk 

From Eqs. (|33|) and (|40)l . we see that we need to evaluate 



E 



6"(k) 



E 



J (fc) 



(53) 



To do so, we note that our phase shifts are those obtained from diagonalizations of our matrices of W 
and W, respectively. In each case, the sum of the tangents of their phase shifts forms a trace which 
is, of course, invariant under changes in the representation. 
Thus we can, for example, write 



h tan(^)] = Trace(W) = £ W j 



(54) 



K 



but since we expand the phase shift in powers of k, and keep at most the first two terms , we find 
that we can evaluate Eq. (|53[) in first order Distorted Born, as: 



h 2 24 y/2 
mg 7r 



E 

.0,3,6.. 



(K + 3) 2 



A(K + 3) s 



- 1 r 

1 3 ^ 



{K + lf 



4(^ + 1)2 - 1 



(55) 



where we have combined the cosine and the sine contributions. 

As we can see, and as is appropriate for our clusters, the individual contributions from the (123) and 
the (12,3) terms are infinite! We expect the difference to be finite. 

Let us comment: 

As shown in a semi-classical treatment [7], each of the phase-shift sums corresponding to the above, 
in 3 dimensions, is associated with a classical expression integrand, which integrated over the position 
variables, diverges as the volume tends to infinity. When all of the terms of the cluster are taken 
together, the resulting integral is finite and volume independent. 

Abel summation 

We see, therefore, that the infinities are the result of taking the infinite volume limit ... and the 
differences of the phase shift sums must be taken with great care. 

The literature of the evaluation and manipulation of divergent series is rich and of great intellectual 
interest. We give some references [20J. 



A famous method to evaluate such series is due to Abel. Using his method we introduce an ad- 
ditional parameter l x\ and powers of x, which will allow us to obtain finite results for each of our 
divergent sums (for < x < 1), subtract our sums, let x = 1 and then revert to our 'physical' results. 



We first start by rewriting the sums in (|5"5|) as 

K 2 



E 



3,6.. 



4 K 2 



1 ^ 
1 3 ^ 



K 2 



4 K 2 - 1 



Following Abel, we multiply the terms of the diverging sums by the strongly decreasing function x K 
where < x < 1 

oo 2 oo 2 

^ 4K 2 ~1 3 ^ 



3,6. 



AK 2 - I' 



1G 



Here we see that for equal values of K, on both sides, say K = 3, the same term appears in b< 
sums, multiplied by x to the same power. This is important. 
Changing the index in the first sum, we obtain 



sr A i k K 

4K 2 - 1/9 3 4 K 2 - 1 



1,2... 

For the first term, then, since the series is absolutely convergent for < x < 1, we can write: 



Kt. xK+3+ kt.. xK+a ((^+3) 2 -i) ( 



Similarly, for the second sum, involving the (12,3) terms, we write: 



- J2 * K+1 + ^ E xK+1 ( ( x + i)2-i) 
0,1,2... 0,1.2... v ^ ' > 4/ 



The expression in the bracket in (|55[) then becomes: 



1 



(^ + 1)2-1 



Collecting the divergent parts, gives us 



0,1,2... v" 1 4 

(• 



1 y ^A'+l _ x " _ 1 g _ X + 2x 

3 f-' 1-z 3 3 1 - a; 3(1 + 2; + a; 2 ) V ' 

0,3,6... 0,1,2... ^ 1 1 / 



which tends to —1/3 when x 1" 
On the other hand 



V - = - V 1 -2 - 

o,fe..^ + 3 ) 2 -I 9 i&.^' 2 -M 7(3) 
and 

v - = v 1 -2 

4-? (k + i) 2 -\ ,4m 7 2 -l " 



0.I..2... (-^ + 1) i.j.:; 

So that we obtain for (15 



m 5 it \ 4 3 16 ^ y(3) 3)1 mg 2 y 
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The oscillatory terms 



Abel summation applied to the sums for the oscillating terms. 



Let us use the Abel procedure for the elements of the sums that appear in the integrands of the 
oscillating terms. I.e., recall that in and (|42p we obtained a term such as (=f)"[cos(2<5 v (A;)) + 1], 
which gives a —2 for the even K terms and a 2 for the odd K terms, as k — > 0. Using, then a 'trace' 
argument, we pass from the sums over v to sums over the indices K associated with our solutions. 
I.e., if we take one of the W - associated with one of the sums in the oscillatory terms - say for a (123), 
even basis - and diagonalize it, we will then obtain S K for a set of values of K, such as 0, 6, 12, etc. 
We can use these values of K to be our values of v, i.e. to be the values of the index over which we 
sum. 

We will then have 



E xK+z - E xK+3 - E xK+1 + E xK+1 

0,6,12... 3,9,15... 0,2,4... 1,3,5... 



which gives 



then 



c 3 (1 — x 3 ) x(l — x) x 3 



1-x 6 1 - x 2 l + x 3 l + x 

x 3 x(l — x + x 2 ) x(x — 1) 
l + x 3 l + x 3 l + x 3 



which goes to as x — > 1. 



Further, if in the integral one makes the change of variable u = kp max , then when one takes the 
limit p m ax -+ oo, the argument of the 5 v (u/ p max ) will go to zero, and if we take the appropriate sums 
we find ourselves with the same expressions that we have evaluated above. The integrand of the total 
oscillatory contribution will be exactly zero. 
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The Fugacity Coefficient 
our result 

Our fundamental formula is: 



(61) 



Putting the values that we have just obtained ((60 j) for the derivative of the difference of the phase 
shift sums, we obtain: 

\ T y nmg 2 J J 
and remembering that At can also be written as (2TTh 2 (3/m) 1 ^ 2 we obtain our result: 



which agrees with the first term of the expansion of the integrals of Dodd and Gibbs[4j. 
See Appendix C. 
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Conclusion 



Our aim, our hope, is to evaluate, express the virial coefficients in terms of scattering, asymptotic 
quantities of the wave functions, and bound state energies. 

This was done in the 1930's for the 2nd virial by Uhlenbeck and Beth[2T], and by Gropper|22|. Our 
dream has been to do the same thing for the 'higher' virials. We believe that we have made progress 
in doing so for all the 'higher' virials - since our approach can be generalized to accommodate any 
number of particles - but demonstrate our method for a prototype third fugacity cluster, required for 
the 3rd virial. 

In this paper, we present results for a model in 1-dim. (with delta function interactions) because 
this model has attracted the attention of many workers[2]l [231 due to its simplicity, which allows 
the evaluation of many results, expressions, and steps, analytically. Further, Dodd and Gibbs[4j, 
uniquely, were able to look at the Statistical Mechanics and evaluate the third fugacity cluster in 
1-dim. and we anticipated that we would be able to compare with their result. 

In essence, our overall reasoning is simple and, we think, transparent! 

Our fugacity clusters are expressed in terms of traces of terms such as Tr(e^ l3Hn ). We know that 
we can introduce complete sets of position eigenstates in these traces, which will then require us to 
integrate over position variables, and a wave number coordinate, as done for our model by Dodd and 
Gibbs|lj. 

We do introduce complete sets of states, which we choose to expand in a hyperspherical adiabatic 
basis, and then integrate analytically over all position coordinates - including our hyper-radius 
coordinate p - obtaining results characterized, here, by the asymptotic behaviour of the wave functions 
(i.e. for large p), described by eigenphase shifts, functions of the wave number. 

This approach is general. Yes, there are subtle problems involving bound states, but here, in our 
1-dim. repulsive case, we need not face them. In 1-dim., however, we find that there is a possibility of 
a contribution from oscillating terms, which for example do contribute[16J to 62, in the case of Bose 
statistics. Our results indicate that there is no such contribution for 63 in Bose statistics. 

A delicate problem that we face, here, is the subtraction between phase shift sums. As we have 
shown, after evaluating the cm. contribution, we must subtract contributions proportional to L 2 , so 
as to obtain L-indcpendent results. In our formalism, the first subtraction takes place, eliminating 
one power of L and yielding 2 types of W-phase shifts. We then saw that we had to subtract these 
two phase shifts sums to yield a finite, i.e. L-independent, answer. Using the results only of order k, 
we were successful in doing so, using a summation procedure associated with Abel. 

For our 3 particles in one dimension, subject to delta- function interactions, we were able, for the 
fully interacting case (123), using a distorted Born approximation, to obtain the phase shifts (or their 
tangents), to order k 2 . In our hyperspherical adiabatic formalism, they are associated with the 'tail' 
of long ranged potentials. 

For the (12,3) case, we encountered reliable results for the lowest contribution, to order k. However, 
for K = 0, and the coefficient of k 2 , we obtain divergences from the contributing terms. Added 
together, these divergences cancel but yield a finite term that we cannot trust. See, in Appendix D, 
our analysis of the domain of applicability of the distorted Born. [We note that we tried - very hard, 
in many ways, using different methods - to obtain the coefficient of the k 2 term, but could not.] 

Using, then, the result of order k, which we do believe is correct in our distorted Born formalism, 
we obtain a result, for the temperature behaviour of 63, which agrees with the leading term in an 
expansion of the double integral proposed by Dodd and Gibbs. 

We feel justified, therefore, in believing in the validity of our formalism, and of its extension to 
1-dim. and Bose particles. 

Ultimately, in our future plans, we will wish to extend, rigorously, our general formalism, to include 
cases of three-body systems with attractive potentials and bound states. Specifically, we need to be 
able to accommodate states in which a possible outcome, asymptotically, would involve, for example, 
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a 2-body bound state + an asymptotically free particle. In an old paper |3], we proposed a formula to 
do this, a very intuitive formula, without a proper proof. Surely, one of our aims must to be remedy 
that! 

(We note that to take into account a 3-body bound state is trivial; in the evaluation of the trace we 
simply add an exponential term involving the bound state energy.) 
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Appendix A 

See the discussion in the article by Larsen and Poll |26| . 

0} (k, p) ~ {kpf 2 {8} Jj (kp) + fl}(fc) N, (hp)) (64) 

As the reactance matrix \\Rj\\ is real and symmetric, it can be diagonalized by a real orthogonal 
matrix Cj- . Consequently, we can introduce eigensolutions: 

# (M ~ (M V2 q (Jj (*p) + A„(fc) A, (M) (65) 

characterized by a unique eigenvalue for all of its amplitudes. 
The eigenphase shifts are then defined by: 

A„(fc) = -tan^(fc) (66) 
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Appendix B 



In this Appendix we give the asymptotic expression of the interaction ( )^ entering the coupled 
equations (Eqs. (|43l44p ). It consists of a series in terms of 1/p and is used in our first and second 
Born approximations. We put for the (123) case 

q K = (K + 3)+J2% (67) 

j 

with K — modulo 6 for K even, and K = 3 modulo 6 for K odd, and for the 2 + 1 case 

qK = {K + l) + Y J '^ (68) 

j 

with K = modulo 2 for K even, and if = 1 modulo 2 for K odd. This leads to (see the previous 
discussions on the calculations of the adiabatic eigenfunctions) 

« - C* + *> - 36 ( | + + 1296 (§±1) + O (69, 
for every value of X in the (123) case and 

^ = v + v - 12 ( + 144 (S?) + ° (?) (70) 

for every value of K in the (12,3) case. 
In both cases c = (2m/h 2 )(3g/TT\/2). 

Obtaining the asymptotic expansion of Ak, C and D matrices (See definitions in Eqs. (|6l9ll0p ) 
in terms of 1/p, we recover (and extend) the results for (123) of Larsen and Popiel |10|.: 

A tf M = W-l/4)/p^ (g + 3 ;- 1/4 +0 (l) (71) 

K > ( _ _d_ lgVg (K + 3)(K' + 3) cos(7r(ir - K')/6) 
k[P) ~ p 2 (K - K')(K + K' + 6) 7T 

<P 1U(K + 3)(K' + 3)cos(n(K ~ K')/6 (\ 



p 3 (K — K'){K + K' + 6) ir 2 \p 4 



° 7T (72) 



and 



K >, _ d_ 2AV2(K + 3)(K' + 3) cos(7r(Jr - K')/6) 

k{P) ~ p 3 (K - K')(K + K' + 6)n 

(f 1U(K + 3){K' + 3)(5(K + 3) 2 - (K' + 3) 2 ) cos^if - K')/6) f JL_ 

p 4 (Jf- A") 2 (A- + A^ + 6) 2 tt 2 + ^p 5 

and for K = K' 

[Note that C|f is an antisymmetric matrix and that, therefore, the diagonal elements are zero.] 
In these expressions d = ti 2 /mg. 
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We obtain in the new case, (12, 3), 

p 2 \p 

K >, , _d_ 4V2 (K + 1)(K' + 1) cos(7r(if - K')/2) 

K [P) ' p 2 [K - K')(K + K< + 2) 7T 

£_ 16(K + l)(K' + l)cos(ir(K-K')/2) 
p 3 {K - K')(K + K> + 2) 7T 2 VP 4 

and 



- _d_ 8V2(JC + 1)(K> + 1) cos(7r(^ - K')/2) 

K [P) ~ p 3 (K - K')(K + K' + 2)tt 

16(K + + l)(5(K + l) 2 - (K' + l) 2 ) cos(7r(X - K')/2) 

p 4 (K - K') 2 {K + K> + 2) 2 tt 2 

for K ^ K', and for K = K' we have 

2 <Vq l^K^j ^2„2\ 
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Appendix C 

Expansion of the integrals of Dodd and Gibbs 

For their concluding equation Eq.(18), in the main part of their paper, Dodd and Gibbs[4j write: 

o i>oo />oo 

h - b °= ° (^3/2 / dg / dte -(l/2fi)(s-+e + st) (e -„-ct _ 1)e -cs (7?) 



(2tt)3^' J Jo 
Denoting, the integral above as X, the integral can be rewritten as 

20 ds dtexp(-(s 2 +t 2 + st))(exp{-\(s + tj) -l)exp(-As), 
Jo Jo 



where we have set A = c^/20, and written again as: 

(2/3) ^ jf" ds e- 2 { e— ^ erfc (*±f) - erfc (£) } 



This then leads to the following expansion: 



I/(2 0) ~ + i + ±+o(- 



Taking the first term, we obtain: 



*3 - *S = 7^) 3/2 (2/3) f-^Ll - - A i (78) 



(2tt)3 v /? ; 20^2^; 4ttV2/3c 

Looking then at the two papers, we determine that (3 of DG equals the (h 2 /2m) that we use, and 
that the 2c, used above, is our g(2m/h 2 ). Substituting then, in their result, gives our formula: 

h - 6? 



4irV2 



2G 



Appendix D 



A little model, testing Distorted Born 

Consider the following differential equation: 

d 2 2a 2 



dr 2 (f + ar) 2 



ip(k,r) = k 2 4>(k,r). (79) 



I.e., a Schrodinger equation with a curious potential which, while finite at the origin, tends to 2/r 2 
at large distances. Thus, it mimics the situation that we have encountered in this work, in which the 
'effective' centrifugal term has changed as a function of distance. 
The exact solution is: 

,., . — a 2 (kr) cos(kr) + (a 2 + k 2 + ak 2 r) sin(Ax) 

i>(k,r) = uTTTL — ^ 80 

k^(l + ar) 

which, for r large, becomes proportional to 

^kr | h/i{kr) - N 3/2 (kr) 

which implies that, introducing a W and a phase shift S, we obtain an exact result: W = —k/a, and 
therefore tan(<5) = k/a. There are no terms in k 2 , or fc 3 for the tangent of S. 

For r large, the potential takes the form 2/r 2 — (4/a)l/r 3 + (6/a 2 )l/r 4 ... and using the — (4/a)l/r 3 
in our first distorted Born yields 

^J^drVk^J 3/2 (kr)Vk^J 3/2 (kr)t^l 7 (81) 

which equals 

- 2 ^ k r dZ (j 3/2 ( Z )) 2 ^-- a (82) 

a Jq z a 

We note that the integrand behaves as z at the origin - hence presents no difficulty there - and so 
would a term in the potential of f/r 4 , which would give 1/z 3 and an integrand which would be finite 
at the origin. A term in the potential of f/r 5 would lead to a divergence. 



If we use our distorted Born formulation to 2nd order, we obtain : 

W = -k/a + (k/a) 2 (2tt/5 - 2tt/5) = -k/a 
where the 1st 2ir/5 comes from the 1/r 4 term in the fst Born, and the 2nd 2ir/5 from the 2nd Born. 

We see that in this little model, when the integrals converge, we obtain the correct answer to order 
k 2 . It is the k 3 term in this model (for W) which plays the role of the k 2 term in our (2+1) case. I.e. 
in these cases [model and (2+1)] the relevant integrals diverge. 

We note that the criterion, given by Mott and Massey for the dominance of the tail in determining 
the low k behaviour, is that £ > (s — 3)/2, where s characterizes the inverse power of r being considered. 
Here, our divergent term is associated with I = 1 and s = 5. The inequality is thus violated and the 
Born (and also our distorted Born) has reached the limit of its application.. 
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